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Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Cole¬ 
man’s ensemble IV-representability conditions. Recently, the topic of pure-state TV-representability 
conditions, also known as generalized Pauli constraints, received increased attention following the 
discovery of a systematic way to derive them for any number of electrons and any finite dimensional¬ 
ity of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of 
the pure-state conditions on the results of reduced density-matrix functional theory calculations. In 
particular, we examine whether the standard minimization of typical 1-RDM functionals under the 
ensemble A-representability conditions violates the pure-state conditions for prototype 3-electron 
systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same 
systems and functionals and compare the correlation energies and optimal occupation numbers with 
those obtained by the enforcement of the ensemble conditions alone. 


I. INTRODUCTION 

In Reduced Density Matrix Functional Theory [l| 
(RDMFT), the total energy of an JV-electron system 
is expressed as a functional of the one-body reduced 
density-matrix (1-RDM) 7(r,r') or its eigenfunctions 
4>jair) and eigenvalues rij a , i.e. the natural orbitals and 
their occupation numbers. To obtain the 1-RDM of the 
ground state, the total energy is typically minimized 
under Coleman’s ensemble IV-representability conditions 
(EC) 0, 

oo 

F, rija = Nfj , 0 < rijc < 1. (1) 

o=l 

These conditions ensure that the 1-RDM corresponds to 
an ensemble of fermionic states. The advantage of the 
EC is their simplicity when compared with the pure-state 
IV-representability conditions (PC). The PC are the nec¬ 
essary and sufficient conditions for the 1-RDM to corre¬ 
spond to a pure state and are also known as generalized 
Pauli conditions [0 0 . 

Although some of the PC were known for a long 
timed, a systematic way to derive all conditions for 
a given number of particles N and finite size M of the 
Hilbert space (being equal to the number of natural spin- 
orbitals) has been recently demonstrated by Klyachko et 
al. d, ;3j. Since the number of natural spin-orbitals needs 
to be finite in any implementation, these conditions can 
be derived for all practical purposes. As the conditions 
depend on both N and M, we will refer to them as the 
N-M conditions in the following. In practice, the num¬ 
ber of PC increases rapidly with the number of electrons 
and natural spin-orbitals. For example, there are four 3-6 
conditions while the number of 3-8 conditions is already 
31. All the PC can be written as inequalities for various 


subsets of the occupation numbers. Geometrically, the 
set of conditions defines a M-dimensional convex poly¬ 
tope. Investigations into the relevance of the conditions 
for physical systems have started recently [7l4l0| focussing 
on the question whether the solutions lie on the bound¬ 
ary of the M-dimensional polytope or inside. In this 
work, we discuss the PC in the context of approximate 
functionals in RDMFT calculations. 

Several approximations of the total energy as an ex¬ 
plicit functional of the natural orbitals and the occu¬ 
pation numbers (or the 1-RDM) have been introduced 
in the last decades (TTM28| and employed in electronic 
structure calculations in combination with the EC. In 
these approximations, the exact kinetic and external po¬ 
tential energy functionals in terms of the 1-RDM are em¬ 
ployed while the electron-electron interaction is approx¬ 
imated. Among these approximations a central position 
is held by the Muller functional [Ill, G 3 where 7 in the 
exchange term of Hartree-Fock theory is replaced by 7 1 / 2 
(power in the operator sense). This approximation was 
shown to overestimate correlation energies substantially 
and several remedies were put forward. For example, 
Goedecker and Umrigar [l3j], proposed a modified form of 
the Muller functional where self-interaction terms (same 
index terms) are explicitly excluded from the double in¬ 
dex summations of the Coulomb and exchange like energy 
terms. 

Gritsenko et al. [ 10 ] introduced a hierarchy of repul¬ 
sive corrections based on the separation of natural or¬ 
bitals into strongly and weakly occupied as well as bond¬ 
ing and antibonding orbitals. The proposed BBCn (n 
= 1,2,3) functionals, and in particular BBC3, reproduce 
quite accurately correlation energies and the dissociation 
of diatomic molecules. A modelling of the BBC3 approx¬ 
imation using a fitting formula was introduced by Rohr 
et al. [l5] leading to the functional known as Automatic 
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Correction 3 (AC3). 

Marques and Lathiotakis [lf| (ML) introduced an em¬ 
pirical functional by replacing the square root depen¬ 
dence on the occupations of the Muller functional with 
a Pade-type formula with parameters fitted to reproduce 
the correlation energies of a test set of molecules. An¬ 
other approximation, called the “Power functional”, tar¬ 
geting the study of periodic systems, was proposed by 
Sharma et al. QAI , by substituting the square-root depen¬ 
dence on the occupations of the Muller functional with 
a power a. A value of a = 0.578 was found to opti¬ 
mize the performance for the functional for both finite 
and extended systems jl9| . This functional was shown to 
reproduce quite successfully the band gaps of semicon¬ 
ductors and insulators including transition metal oxides 

Ena. 

Based on the reconstruction of the two-body reduced 
density matrix in terms of the 1-RDM, Piris and cowork¬ 
ers developed a series of approximations [20l - [23| termed 
PNOFn, n = 1,...6 which were proven to be accu¬ 
rate in reproducing several properties of molecular sys¬ 
tems. Functionals that one could classify as extensions of 
RDMFT are those that depend explicitly on the phases 
of natural orbitals [24} like for example the exact func¬ 
tional for 2-electron systems of Lowdin and Shull [30] and 
the theory of the anti-symmetrized product of strongly 
orthogonal geminals (APSG) [25}. APSG accurately pre¬ 
dicts structural and vibrational properties of molecules as 
well as single-bond breaking [26H28j . 

Recently, and in order to address the large numeri¬ 
cal cost of the orbital minimization in RDMFT, local- 
RDMFT was introduced [HI, 32j] which combines the 
non-idempotent 1-RDMs of RDMFT with a local, Kohn- 
Sham-like potential in a variational manner. 

One could argue that, in the minimization of the total 
energy in RDMFT, PC should be enforced in addition 
to the EC. For the exact functional, however, due to the 
minimization principle, the EC are sufficient to guarantee 
that the minimizing 1-RDM corresponds to a pure state, 
i.e. the ground state, in all cases without ground-state 
degeneracy. For approximate functionals, however, this 
is not the case. In other words, in the minimization of an 
approximate functional of the 1-RDM, if one applies the 
EC alone or the PC in addition, there is no guarantee 
that the optimal 1-RDMs will coincide. 

In this paper, we assess how well different approxi¬ 
mate RDMFT functionals, combined with only the EC, 
satisfy pure-state conditions for different 3-electron sys¬ 
tems, namely the lithium atom, the LiH + and the HeJ 
ions, and the hydrogen trimer, using 6 or 8 natural 
spin-orbitals in the calculation. In addition, we perform 
RDMFT calculations for these 3-electron systems enforc¬ 
ing the PC as additional constraints and compare the re¬ 
sults with and without enforcing them. The goal of this 
work is to get a first feeling on whether imposing the PC 
has any impact on the RDMFT results or if the EC are 
sufficient. Three electron systems are the smallest sys¬ 
tems where one can investigate this question because for 


closed-shell systems with time-reversal symmetry the EC 
are proven to be necessary and sufficient for pure-state 
iV-representability [ 6 j]. 

In section hd we briefly review the known conditions 
for pure-state IV-representability for 3-electron systems 
and for Hilbert spaces of 6 and 8 natural spin-orbitals. 
We discuss our results for the four different systems with 
and without enforcing the PC in section Hill and conclude 
our findings in section HVl 

II. PURE STATE IV-REPRESENTABILITY 
CONDITIONS 

Similarly to the EC, the conditions for the pure-state 
IV-representability concern only the occupation numbers. 
To express the PC, at least in the form known so far, 
occupation number ordering according to their value is 
essential. Thus, in order to write the PC we have to map 
the set of occupation numbers {rij a } onto an ordered set 
{A*,}, Ai > A 2 > ... > A m regardless of the spin they 
correspond to. As a realistic example considered in this 
work, for the lithium atom, the Muller functional with 
the EC alone and 6 natural spin-orbitals, we obtain 

m t = 1.000000, n 2t = 0.865779, n 3t = 0.134220, 
n u = 0.999969, n 24 . = 0.000027, n 34 . = 0.000005, 

( 2 ) 

which are then mapped according to 

nit -A Ai, n 2 | -A A 3 , n 3 -|- -a A 4 
n.14, -A A 2 , n 2 ! -A A 5 , n 3 ! -A A 3 . 

We note that using three occupation numbers for each 
spin channel is a question of choice. One could choose 
an unbalanced distribution of the 6 occupation numbers 
into the two channels. For systems of 3 electrons, and 
for a Hilbert space of 6 natural spin-orbitals, there are 
4 conditions which are given in Table |U The first three 
of the conditions are originally inequalities, like all Kly- 
achko’s conditions, with the left-hand-side being less or 
equal to one, however, they reduce to equalities due to 
the ensemble IV-representability condition Eq. 0 which 
the occupation numbers have to satisfy as well. As one 
can see, for the example above, only the first and the 
third conditions are well satisfied while the second is vi¬ 
olated slightly and the forth quite severely. 

The conditions in Tabic U were originally formulated 
by Borland and Dennis in the 70’s [§}. For a 3 elec¬ 
tron wavefunction expressed in terms of 6 natural spin- 
orbitals, 6 , (three for each spin channel) one can 

form 9 Slater determinants that could contribute to the 
wavefunction. Borland and Dennis found numerically 
that only four of these determinants have non-zero co¬ 
efficients. These determinants are |</>i</> 2 <(> 3 |, |</>i04</>5|, 
102 ^ 406 1 and 1 03 05 061 - We label the expansion coeffi¬ 
cients of these determinants with ci, c 2 , c 3 , and C4, re¬ 
spectively. Forming the corresponding 1-RDM we can 
associate every occupation number A$ with two expan¬ 
sion coefficents Ci because each orbital is present in only 
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# 

Condition 

#l 

Condition 

1 

Ai + A6 = 1, 

1 

Ai T A2 T A4 A7 < 2 

2 

A2 + As = 1 , 

2 

Ai + A2 + A5 + A6 < 2 

3 

A3 + A4 — 1, 

3 

A2 + A3 + A4 + A5 <2 

4 

As + A6 — A4 > 0 . 

4 

Ai + A3 + A4 + Ae < 2 


TABLE I. 3-6 PC, i.e. for 3 electrons in a space of 6 natural 
spin-orbitals Q. 
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FIG. 1. Violation in electrons (e) of the 3-6 PC for Li, H 3 , 
LiH + , and HeJ, with different RDMFT functionals minimized 
under the EC alone. On the horizontal axis, the conditions are 
numbered in the order shown in Table [I] A Missing symbol 
implies that the particular condition is not violated by the 
corresponding functional. 
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Ai + X 2 — A 3 < 
A2 + As — A7 < 
Ai + Ag — A 7 < 
A2 + A4 — A6 < 
Aj + A 4 — A 5 < 
A3 + A4 — A7 < 


Ai + As < 1 


A 2 — A 3 — A 6 — A 7 < 0 
A 4 — A 5 — A 6 — A 7 < 0 
Ai — A 3 — A 5 — A 7 < 0 


A 2 + A 3 + 2 A 4 — A 5 — A 7 + As < 2 
Ai + A 3 + 2 A 4 — A 5 — Ag + As < 2 
Ai + 2 A 2 — A 3 + A 4 — A 5 T As < 2 
Ai + 2 A 2 — A 3 + A 5 — Ag -j- As < 2 


Ai + A2 — 2A3 — A4 — A5 < 0 
Ai — A2 — A3 + Ae — 2A7 < 0 


Ai — A3 — A4 — A5 + As < 0 
Ai — A2 — A3 — A7 + As < 0 


2Ai — A2 + A4 — 2A5 — Ag + As < 1 
A3 + 2A4 — 2A5 — A6 — A7 + As < 1 
2Ai — A2 — A4 + Ag — 2A7 + As < 1 
2Ai + A2 — 2A3 — A4 — Ag + As < 1 
Ai + 2A2 — 2A3 — As — A6 + As < 1 


2Ai — 2A2 — A3 — A4 + As — 3A7 + As < 0 
—Ai + A3 + 2A4 — 3A5 — 2A6 — A7 + As < 0 
2Ai + A2 — 3A3 — 2A4 — As — A6 + As < 0 
Ai + 2A2 — 3A3 — A4 — 2A5 — A6 + As < 0 


two Slater determinants. For example, Ai = |ci | 2 + |c 2 1 2 
and \ e = | C 3 1 2 -f- | C 4 1 2 . Since the wave function has to 
be normalized, Ai + Ag has to be one, which is the first 
Borland - Dennis condition. In the same way, one ob¬ 
tains the second and third conditions. Moreover, every 
\ci\ 2 has to be non-negative from which the forth condi¬ 
tion follows. In this way, Borland and Dennis proved the 
sufficiency of their conditions and referred to a work by 
M.B. Ruskai and R.L. Kingsley for their necessity. This 
work, however, was only published 35 years later by M.B. 
Ruskai [ddf. 

The number of conditions increases dramatically with 
the number of natural spin-orbitals that one includes in 
the calculation. Using invariant theory and represenation 
theory, Klyachko and coworkers derived and proved an 
algorithm for computing all such Pauli-like constraints 
for a given number of electrons and natural spin-orbitals. 
For the case of 8 orbitals they obtained 31 conditions that 
need to be enforced which are shown in Table mi In¬ 


in. RESULTS 

We perform all the calculations using the cc-PVTZ 
Gaussian basis set with the HIPPO computer code [35j- 
For the 3-electron doublet state we used the assumption 


TABLE II. 3-8 PC, i.e. for 3 electrons in a space of 8 nat¬ 
ural orbitals. Horizontal lines distinguish different types of 
conditions [ 34 } • 


that the natural spin-orbitals appear in pairs with iden¬ 
tical spatial dependence [36!], i.e. the spin dependence 
enters only through the occupation numbers. This re¬ 
striction is not expected to influence the findings of the 
present work since both PC and EC concern the occupa¬ 
tion numbers only. We performed two kinds of calcula¬ 
tions: (i) standard RDMFT calculations using the Muller 
EM3, BBC3 [3, ML [3 and PNOFO [ij, HI func¬ 
tionals to examine whether, and to what extend, PC are 
satisfied by the optimal occupations, and (ii) RDMFT 
calculations with the 3-6 or 3-8 PC enforced in addition. 
We note that only 3 or 4 spatial orbitals ( 6 , 8 natural 
spin-orbitals) were assumed to have non-zero occupation 
even in case (i) in order to conform with the Hilbert space 
restrictions of the PC and to allow for fair comparisons. 
We also note that the large size of the basis set allows 
for the use of the optimal M-dimensional Hilbert space 
in each calculation. Occupation number and orbital op¬ 
timizations were performed with customized sequential 
quadratic and conjugate gradient methods, respectively. 
Our calculations are for testing purposes, since 6 spin- 
orbitals are generally not enough to describe dynamic 
correlations in 3-electron systems. However, they are still 
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reasonable because for 8 spin-orbitals the smallest occu¬ 
pation number is typically of the order of 10 ~ 5 - 10 ~ 6 . 

We first examine whether the occupation numbers ob¬ 
tained with standard RDMFT calculations, i.e. by en¬ 
forcing only the EC, satisfy the PC. In Fig. [T] we show 
the violation of the 3-6 PC in electrons for Li, LiH + , 
HeJ, and H 3 . As the first 3 of the 3-6 conditions are 
equalities the quantity we refer to as “violation” is the 
sum of the relevant occupation numbers minus 1 (the rhs. 
of all equalities). The fourth condition is violated when 
A 5 + A 6 — A 4 < 0 and the violation is the difference of the 
lhs. of the condition from zero, i.e. the minimum value 
of the lhs. satisfying the condition. In a similar way, we 
determine the violation of the 3-8 conditions. For the 
3-6 conditions, we found that, for all functionals and all 
systems apart from H 3 , the fourth condition is always 
violated and has the highest violation. While for H 3 all 
violated conditions are of the same order of magnitude, 
for the rest of the systems the fourth condition is vio¬ 
lated the most by at least an order of magnitude. In 
all cases, the ML functional appears to give the smallest 
violations among the approximations employed. Inter¬ 
estingly, when both the first and the second conditions 
are violated, as in Li and LiH + , they are violated by 
the same amount. Due to the enforcement of the en¬ 
semble IV-representability conditions and the pinning of 
the largest occupation number to one, the mapping ([3j) 
implies that both the first and the second conditions re¬ 
duce to n. 34 , = 0. Hence, they are violated by the same 
amount. For the hydrogen trimer and the He((~ molecule, 
the mapping onto ordered occupation numbers is differ¬ 
ent than in Eq. ©■ As a result, the first and the third 
conditions reduce to 7134 . = 0 and are, hence, equally vi¬ 
olated. For Li and LiH + , due to the pining of the first 
occupation number, Ai = 1 , the third condition is au¬ 
tomatically satisfied since the total number of electrons 
with spin up equals Ai + A 3 + A 4 = 2. 

In Fig. [2j we show the violation of the 3-8 conditions. 
As we see, the conditions with numbers 9, 10, 20, and 
28 in Table Qj] are always satisfied while those with num¬ 
bers 5, 18, and 26 are always violated. For weakly cor¬ 
related systems, where Ai to A 3 are close to one and all 
remaining occupations are close to zero, it is clear that 
conditions 20 and 28 will always be satisfied. Condition 
5 is always violated because we find that Ai is pinned 
to one for all systems and functionals considered here, 
which implies that the only way to satisfy the condition 
is A 2 = A 3 which is actually not the case. For the set of 
3-8 conditions, again the ML functional gives the small¬ 
est violation, with PNOFO coming next, while between 
the Muller and BBC3 functionals there is no clear order. 
These two functionals, however, show a similar violation 
for many of the conditions for all systems except H 3 . Due 
to the larger number of conditions there are also several 
conditions that are violated to a similar extend, i.e. con¬ 
trary to the 3-6 case there is not one dominant violation. 

In order to assess the impact of the PC on RDMFT 
calculations, we performed energy minimizations imple- 


A. For the 3-6 conditions: 



Muller 
\E^\ A E c 

BBC3 
|M B) | A Ec 

ML 

|^ E) | A E c 

PNOFO 
|^ E) | A Ec 

Li 

12.2 

11.0 

10.0 

9.2 

2.4 

1.2 

2.8 

2.0 

h 3 

45.5 

2.1 

24.1 

2.5 

17.1 

0.5 

20.2 

0.5 

LiH+ 

16.1 

1.5 

7.7 

5.6 

19.1 

15.5 

11.8 

0.5 

He+ 

81.7 

39.1 

71.4 

43.3 

25.8 

2.0 

32.8 

7.5 


B. For the 3-8 conditions: 



Muller 
|4 E) | A Ec 

BBC3 
|fi£ B) | A Ec 

ML 

\Ei E) \ A E c 

PNOFO 
|Fc E) | a e c 

Li 

23.4 

20.0 

16.5 

15.2 

4.8 

1.3 

5.6 

3.6 

h 3 

53.8 

0.7 

29.1 

2.3 

23.3 

0.2 

25.0 

0.4 

LiH+ 

19.9 

1.0 

10.3 

6.2 

22.2 

1.8 

15.2 

3.2 

He+ 

90.7 

38.0 

76.6 

42.4 

34.3 

1.3 

37.9 

6.3 


TABLE III. Absolute correlation energy, |i?c E ' > |, for the EC 
using 6 or 8 natural spin-orbitals, and its change, A E c , when 
the 3-6 or 3-8 PC are enforced in addition. All energies are 
given in mHa. 


menting the PC as additional constraints. A sensible 
comparison is that of the correlation energies Ec ; for 
the EC alone and Ec when the PC are additionally 
enforced, defined as = E — i?( HF ) and E^ = 
£( p ) - £( hf \ where E^ e \ E^ are the total energies 
obtained with the enforcement of the EC alone, and the 
PC in addition, respectively, and £A hf ) is the Hartree- 
Fock total energy. The enforcement of the PC as addi¬ 
tional constraints raises the total energy by the difference 
A E c = Ep -E^ which will be zero only in those cases 
where the minimization without PC enforcement fulfills 
all the PC. No such case was found, however, in the sys¬ 
tems and functionals we considered. In Table mu we 
include the absolute correlation energies |I?c | without 
the PC, as well as A E c for the different systems and func¬ 
tionals. Enforcing the PC for the Li atom, and for HeJ 
raises the correlation energy by an amount A E c which is 
of the same order of magnitude as the correlation energy 
E& itself, for almost all functionals. For the other two 
systems, the effect of the PC on the correlation energy 
varies depending on the system and the functional but 
it is generally smaller. In general, however, A E c is quite 
substantial. As expected, the ML functional which was 
found to violate the PC less than the other functionals 
yields smaller A E c . The trends we observe are similar 
for the 3-6 and the 3-8 PC. 

State pinning, i.e. the existence of natural spin-orbitals 
with border occupations (1 or 0 ) is very common within 
RDMFT approximations. More specifically, with the en¬ 
forcement of the EC alone there exist many spin-orbitals 
(core states) that are pinned at one while pinning at zero 
is rather unlikely, i.e. there is an “over-pinning” of core 
states at one for many approximations. However, this 
inaccuracy is not a severe flaw in most cases, since these 
natural spin-orbitals have an occupation very close to one 
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FIG. 2. Violation in electrons (e) of the 3-8 PC for the Li atom, and the H 3 , LiH + , and He^ molecules with different RDMFT 
functionals minimized under the EC alone. On the horizontal axis, the conditions are numbered in the order shown in Table HT1 
A Missing symbol implies that the particular condition is not violated by the corresponding functional. 


A. For the 3-6 conditions: 



Muller 

BBC3 

ML 

PNOFO 

Li 

9.6-10 -3 

1.7-10 -7 

l.MCT 4 

3.6-10 -5 

h 3 

4.2-10 " 3 

4.8-10 -4 

8.5-10 -4 

1 . 1 - 10 -3 

LiH+ 

1.3-10 -4 

0 

3.4-10 -5 

1.1T0 -6 

He+ 

4.1-10 -4 

9.6-HT 5 

3.4-1CT 4 

1 . 6 - 10 -4 


B. For the 3-8 conditions: 



Muller 

BBC3 

ML 

PNOFO 

Li 

6 . 1 - 10 -4 

3.0-10 -4 

1.5-10 -4 

2 . 2 - 10 -4 

h 3 

1 . 1 - 10 -3 

2.9-10 -4 

5.0-10 -4 

3.6-10 -4 

LiH+ 

9.1-10 -4 

4.6-10 -4 

7.8-10 -4 

9.4-10 -4 

He+ 

3.610 -3 

1 . 6 - 10 -3 

6 . 0 - 10 -4 

1.5-10 -3 


TABLE IV. Deviation from pinning for the largest occupa¬ 
tion number, 1 — Ai, with the 3-6 or 3-8 PC enforced. With 
the enforcement of the EC alone, this occupation number is 
pinned to 1 in all cases considered. 


anyway, and they are usually core states not affecting 
the chemical bonding. Although the existence of pinned 


states in the exact 1-RDM of interacting electron sys¬ 
tems is not an easy question, simple 3-electron systems 
like those studied in this work should not have pinned 
states. However, with only the EC enforced, for all the 
functionals employed here Ai = 1 for all systems. Inter¬ 
estingly, as we see in Table m applying the 3-6 PC in 
addition, Ai becomes fractional in all cases apart from 
one (LiH + with the BBC3 functional). Moreover, im¬ 
posing the 3-8 conditions leaves no pinned states and in 
most cases unpins Ai even more than the 3-6 conditions 
for the corresponding system and functional. This is due 
to enforcing condition 5 which, as discussed above, can¬ 
not be satisfied if Ai = 1 unless A 2 = A 3 . As the latter is 
not the case in any of our systems the only way to satisfy 
this condition is the unpinning of Ai. 

In addition to artificially pinning states, many approx¬ 
imate functionals are also known to overestimate the oc¬ 
cupation of the weakly occupied states [37]. Due to the 
unpinning of the largest occupation number, one might 
expect that the occupation of those states is further in¬ 
creased. However, for both the 3-6 conditions and the 3-8 
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Muller 

BBC3 

ML 

PNOFO 

Li 

1, 5 

2, 5, 6, 16 

5 

1, 5 

h 3 

2, 5, 7, 11, 
18, 25, 26 

2, 5, 7, 11, 
18, 25, 26 

5, 11, 18, 26 

4, 7, 11, 18, 
23, 25 

LiH+ 

1, 2, 4, 5, 
7, 18 

1, 2, 3, 4, 

5, 6, 7, 8, 

10, 12, 14, 16, 
17, 18, 31 

2 , 6, 10 

2, 3, 6, 7 

hST 

2, 5, 7, 18 

2, 5, 7, 18 

2, 5, 7, 18 

2, 5, 7, 18 


TABLE VI. Saturated 3-8 conditions, i.e. those conditions 
where the set of the optimal occupations lies on the boundary 
of the polytope when enforced. For the numbering convention 
see Table [III 


A. For the 3-6 conditions: 



Muller 
w/o with 

BBC3 
w/o with 

ML 

w/o with 

PNOFO 
w/o with 

Li 

13.42 

0.11 

11.64 

0.09 

0.17 

0.08 

0.61 

0.04 

h 3 

10.89 

9.85 

6.41 

5.26 

0.52 

0.49 

1.93 

1.91 

LiH+ 

1.04 

0.33 

1.09 

0.15 

0.18 

0.12 

0.42 

0.17 

He+ 

14.56 

3.72 

13.78 

2.86 

0.54 

0.35 

2.26 

1.20 


B. For the 3-8 conditions: 



Muller 
w/o with 

BBC3 
w/o with 

ML 

w/o with 

PNOFO 
w/o with 

Li 

22.23 

0.29 

15.53 

0.12 

0.35 

0.18 

1.36 

0.10 

h 3 

11.79 

11.30 

6.50 

5.33 

0.63 

0.62 

2.22 

2.21 

LiH+ 

1.51 

0.67 

1.24 

0.18 

0.27 

0.25 

0.53 

0.33 

He+ 

15.03 

4.23 

13.85 

2.79 

0.64 

0.46 

2.34 

1.33 


TABLE V. Sum of the occupation numbers (xlOO) of weakly 
occupied orbitals with or without enforcing the 3-6 or 3-8 PC. 

I 


Since all the 3-8 constraints are inequalities one distin¬ 
guishes between those constraints that lie on the bound¬ 
ary of the polytope when being enforced, i.e. where the 
lhs. is equal to the rhs. and those constraints that lie 
inside the polytope, i.e. where the lhs. is truly smaller 
than the rhs. Physically, being on the boundary implies 
that any further reduction of the total energy will come 
at the cost of violating at least one of these constraints. 
In Table IVll we list those conditions that are lying on the 
boundary for each system and functional. As one can see, 
the smallest number of boundary constraints appears for 
the ML functional which also showed the smallest viola¬ 
tion of the constraints when they were not enforced. 


IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we examined if the additional enforce¬ 
ment of the pure-state 7V-representability conditions, also 
known as generalized Pauli constraints, has an impact 
on RDMFT calculations. More specifically, we checked 
if the standard RDMFT optimization for a few repre¬ 


conditions (see Table [V]) this is clearly not the case. For 
all systems, and all functionals, the occupation of the 
weakly occupied orbitals is actually decreased in some 
cases even quite significantly. 


sentative approximations under the EC of Coleman re¬ 
spect the PC conditions for 3-electron systems, namely 
the Li atom, the linear Hydrogen trimer, and the posi¬ 
tive ions of the He 2 and LiH molecules. We found that 
in all cases at least some of the pure-state conditions are 
violated. We then applied the 3-6 and 3-8 PC as addi¬ 
tional constraints in the minimization and found that the 
increase in the total energy corresponds to a significant 
fraction of the correlation energy obtained with the EC 
alone. In addition, we found that the optimal 1-RDM 
is qualitatively improved by the enforcement of PC since 
the pathology of state pinning is alleviated and the to¬ 
tal charge occupying weakly occupied orbitals is reduced. 
Finally, we found that, in analogy to state pinning of the 
EC, many PC are saturated, i.e. the optimal 1-RDM lies 
on several boundaries of the polytope of pure-state N- 
representability. Thus, in conclusion, we can claim that 
the enforcement of PC has a significant impact on the re¬ 
sults of RDMFT approximations that could potentially 
be explored to define better approximations. 

Unfortunately, the number of N-M PC for large num¬ 
ber of electrons N and sizes of the Hilbert space M can 
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be quite large and difficult to handle in practical imple¬ 
mentations. However, the existence of a systematic way 
to derive these conditions for arbitrary N and M and 
their linear form should facilitate their adoption in more 
general applications of RDMFT in the future. 
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